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A generalization of the operator product expansion for Euclidean correlators of 
gauge invariant QCD currents is presented. Each contribution to the modified 
expansion, which is based on a delocalized multipole expansion of a perturbatively 
determined coefficient function, sums up an infinite series of local operators. On 
a more formal level the delocalized operator expansion corresponds to an optimal 
choice of basis sets in the dual spaces which are associated with the interplay 
fH , of perturbative and nonperturbative Af-point correlations in a distorted vacuum. 

^ ^ A consequence of the delocalized expansion is the running of condensates with 

I ' the external momentum. Phenomenological evidence is gathered that the gluon 

condensate, often being the leading nonperturbative parameter in the OPE, is 
indeed a function of resolution. Within a model calculation of the nonperturbative 
' corrections to the ground state energy of a heavy quarkonium system it is shown 

exemplarily that the convergence properties are better than those of the OPE. 
Potential applications of the delocalized operator expansion in view of estimates 
. of the violation of local quark-hadron duality arc discussed. 

M. 

1. Introduction 

1.1. Some historical remarks 

Once again: Happy Birthday, Arkady! This is a very hvely conference that 
I am happy to attend. 

In this talk I would like to report about work done jointly with Andre 
Hoang. The subject is to formulate a framework for a generalization of Wil- 
son's operator product expansion (OPE), to confront this framework with 
experiment, and to demonstrate its usefulness in some model calculations. 

As a warm-up I remind you of some of the OPE applications. The OPE 
is a powerful tool for the treatment of conformally invariant field theories 
where its convergence is explicitly shown. Unfortunately, this feature does 
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not survive in realistic, four-dimensional field theories. Here the expansion 
is believed to be asymptotic at best. 

The construction of effective theories like HQET or of the effective 
Hamiltonians of electroweak theory are crucially relying on the OPE which 
determines effective vertices by a perturbative matching of amplitudes with 
those of the fundamental theory. Heavy fields are effectively integrated out 
in this procedure (sec [ 1] and references therein). 

The OPE turned out to be important for a semi-perturbative derivation 
of the exact (3 function m N =1 super Yang-Mills theory [ 2] . 

And, finally, there is a major application in QCD sum rules [ 3] where 
the "theoretical side" of a dispersion relation for the correlator of gauge in- 
variant currents is represented in terms of an OPE. The occurence of power 
corrections to the low-order in perturbative result can be argued from 
the summation of bubble graphs inducing so-called infrared renormalons, 
and formal contact to the OPE can be made [ 4] . A practical version of 
the OPE was proposed by Shifman, Vainshtein, and Zakharov in [ 3] where 
nonperturbative factors in the power corrections, the QCD condensates, 
are viewed as universal parameters of the QCD vacuum to be determined 
by experiment. Due to the universality of the condensates this method is 
rather predictive in practice. 

1.2. Problems of the OPE in QCD 

Despite its many successes in the framework of the QCD simi rule method 
there are some unanswered questions concerning the convergence properties 
of the OPE. Answering these questions very likely resolves the problems 
connected with local quark-hadron duality (see [ 5] in references therein) 
which can be formulated as follows. 

Given analyticity of the correlator away from the positive real axis and 
invoking the optical theorem the discontinuity of the correlator at time- 
like momenta should yield the spectral function of the hadronic channel 
which corresponds to the QCD current in the correlator. This procedure, 
however, fails badly if a "practical" OPE with a small number of power 
corrections is used. Whereas the experimentally measured spectral function 
shows resonance wiggles at moderate time-like momenta a continuation 
of the OPE generates a smooth behavior down to small momenta. This 
may not harm the sum rule method since the relevant quantity there is a 
spectral integral, that is, an average over resonances which agrees quite well 
with the average over the associated discontinuity of the OPE. However, in 
applications, where the local properties of OPE discontinuities are needed 
[ 6] , the violation of local quark-hadron duality can induce large errors. 
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Another concern, which is connected with the violation of local quark- 
hadron duality, is the asymptotic nature of the expansion. Asymptotic 
behavior can be argued from estimates of operator averages using the in- 
stanton gas approximation [ 3] or by appealing to the renormalon idea. 
From the OPE itself, where naively an expansion in powers of Aqcd/Q 
is assumed, there is no way of estimating the critical dimension at which 
the expansion ceases to approximate. In the framework of a dclocalizcd 
expansion it was argued in [ 7] that the knowledge about nonperturbative, 
gauge invariant correlation functions would allow for such an estimate if 
the decay of effective correlation length with mass dimension is sufficiently 
fast. 

The central theme of this talk is a nonlocal generalization of the OPE 
which has the potential to cure the above short- comings. 

2. Delocalized operator expansion 

The OPE of the correlator of a current j reads 

i j d^xe'*- {Tj{x)m) = Q^^X^CdKQ') {Odim , iQ' = -q' > 0) , 

d=0 1=1 

(1) 

where d, I run over the dimension and the field content of the local operator 
Odi, respectively. The Wilson coefficient Cdi is perturbatively calculable. 
The expansion (1) contains reducible chains of operators, such as 

GD'^G, GD'^G, GD^G... , 
qq, qDi,^fj,q, qD'^q, qD'^Df^-y^.q , (2) 

Dn denotes the gauge covariant derivative. Here "reducible" refers to the 
fact that each chain can be obtained from a Taylor expansion of an as- 
sociated gauge invariant 2-point correlator which is not decomposible into 
gauge invariant factors. Examples, which correspond to the chains in (2), 
are the gluonic field strength correlator [ 8] 

and the bilocal quark "condensate" 

where P denotes path ordering. For lattice measurements the path connect- 
ing the points and x is taken to be a straight line. This should capture 
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the largest scale A in the nonperturbative decay of the correlator. The gen- 
eral case of A/'-point correlations will be addressed below. The effect of a 
delocalized expansion is a partial summation of such a chain in each order 
of the momentum expansion of the perturbative coefficient. This leads to 
an improvement of the convergence properties. 

Let me make this explicit for the case of 2-point correlations in a 1- 
dimcnsional Euclidean world [9]. The contribution to the OPE, which 
arises from this, can formally be written as 



/ dxf{x)g{x)= dxdyf{y)d{x-y)g{a 

J —OO J 



= j dxdyfiy) 

= E [ / dyf{y)y''] \f dx t})ls^^^ {x) g{x 



■OC' 

E 

n=0 



fn{oo)g„{oo) . 



(5) 



Here f{x) is associated with perturbative correlations whereas g{x) de- 
scribes nonlocal, nonperturbative effects. The integral expressions in the 
first and second square bracket correspond to local Wilson coefficients and 
operators of the like of (2), respectively. From (5) it is seen that the part 
of the OPE being generated by 2-point correlations is induced by a bilinear 
form (•,•) = J dx ■ X ■. The expansion (5) corresponds to a specific choice 
of basis {e„(oo), e„(oo)} in the associated dual space of square integrable 
functions where 

(-1) 



e„(oo) = y" 
and we have orthonormality 



and e„(oo) 



n! 



■S^"\x). 



( ; ) 



From Eq. (5) it follows that 



E 

n=0 



/ 



dy f{y) y" 



i-ir 



(6) 



(7) 



(8) 



This corresponds to an expansion of the perturbative factor f{x) into de- 
generate, that is, zero-width multipoles. It is suggestive that an expansion 
into multipoles of width ^ which is comparable with ~ Q~^, con- 
verges better (see Fig. 1). Like the OPE this expansion is controlled by 
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x=0 X 
Figure 1. Schematic drawing of the short-distance function f(x) and the nonperturba- 
tive 2-point-correlator g{x) illustrating the scale hierarchy Ay/Ag ~ A/Q ^ 1. 



powers of ^ ~ ^, and it breaks down if A/ ~ Deltag. The expansion 
into finite-width multipoles corresponds to a change of basis in dual space. 
Rather than expanding into the basis {e„(oo), e„(oo)} one may expand into 
{e„(0),e„(0)} where 

et(^) = ^i|#, e^W-^i^n(0.)e-^^^^ (9) 



We then have 



/ 



dx fix) g{x) = fni^) 9n{^) ■ (10) 

n=0 



In (9) Hn denotes the nth Hermite polynomial. An expansion into any 
other basis that allows for a sensible multipole expansion of f{x) is possible, 
and the restriction to a Gaussian based expansion is just for dcfinitcncss. 
Note that (9) reduces to (6) in the limit Q ^ oo. According to (9) the 
transformations, which link a basis with resolution parameter Q to a basis 
with f2', form a 1-parameter group. The moments fnif^) in (10) have finite 
expansions into the Wilson coefficients fn{oo) 

/0(O) = /o(oo) , /2(0) = /2(00) - ^/0(00) , • • • . (11) 

Chosing a particular basis and a value for Cl is analogous to a choice of 
regularization and normalization point in pure perturbation theory. In 
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practice (11) enables to consider perturbative renormalization and delo- 
calized operator expansion (DOE) separately. The expansion of g{Q,) into 
local condensates is inSnite, for example 

50 (fi) =50(00) + ^52(00) +••• . (12) 

To evaluate sums like in (12) one would have to know infinitely many 
anomalous dimensions and condensates. On the other hand, one may be 
more pragmatic and start with a (lattice inspired) model for g{x) which is 
believed to contain the information about all (7„(oo). To each order ^ one 
may then derive evolution equations for /„(ri) and 5„(r2) which sum these 
contributions to all orders [ 9] . 

Setting = Q, the parametric organization of the expansion in ^ up 
to an overall factor is 

/o(Q)x5o(Q)~(^)°x{(^)%(^)%-} 
MQ)x..(Q)^(^)\{(^)V(^)%...} 

! . (13) 

The generalization to /C-dimensional Euclidean space is straightforward [ 
9] . It also allows for a treatment of A''-point correlations in 4-dimcnsional 
space. In particular, the n = term of the contribution to the DOE from 
2-gluon correlations yields the local Wilson coefficient times a "runnning 
gluon condensate" {Q = 0,). 

3. Applications 

3.1. Gluonic field strength correlator 

To have a wcll-motivatcd model for g{x) for the case of 2-gluon correlations 
we appeal to a parametrization of the gluonic field strength correlator (3) 
as it was proposed in [ 8] 

In a lattice measurement [ 10] the scalar functions D and Di were 
parametrized as 

D{x'^) = Aoexp[-|a;|/AA] + exp[-|a;|/Aa] , 
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£>i(a;') = Ai exp[-|a;|/AA] + ^ exp[-|x|/AJ , (15) 

Fl 

where the purely exponentially decaying terms are attributed to purely 

nonperturbativc dynamics. The measurements at various quark masses for 
dynamical light quarks indicate that for |x| < the nonperturbative con- 
tribution to D dominates the nonperturbative contribution to Di implying 
that the tensor structure of g^uK\ is identical to that of the correspond- 
ing local operator obtained by letting x ^ 0. This feature leads to great 
simplifications. 



3.2. Extraction of the running gluon condensate 

Using the DOE instead of the OPE, the running gluon condensate can 
be extracted from experiment in channels where 2-gluon correlations are 
dominant. We have done this for the V + A correlator and the charmonium 
system. For a more detailed presentation see [ 9] . 



3.2.1. Sum rules for the V+A correlator 

The spectral function for light quark production in the V + A channel has 

been rcmcasured recently from hadronic r decays for < rn^ by Alcph [ 11 
] and Opal [ 12]. In the OPE the associated current correlator is dominated 
by the gluon condensate, and the dimension 6 power corrections that are 
not due to a double covariant derivative in the local gluon condensate are 
suppressed. The corresponding currents are jf^^^ = U'y^{l ± 75)^ and the 
relevant correlator in the chiral limit reads 

i J d^xe^i- {Tjl:{x)j^{Q)) = {q^q. - q^g^.) n^+^lQ^) , = -q^ . 

(16) 

Since the correlator n(Q^)^+^ is cutoff-dependent itself, we investigate the 
Adler function 



D(Q-) = -Q- dU^^HQ') =91 Tds (17) 



For the corresponding experimental V-l-A spectral function we have used 
the Aleph measurement in the resonance region up to 2.2 GeV^. For the 
continuum region above 2.2 GeV^ we used 3- loop perturbation theory for 
ag{Mz) = 0.118, and we have set the renormalization scale /i to Q. We 
remark that the pion pole of the axial vector contribution has to be taken 
into account in order to yield a consistent description in terms of the OPE 
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Q (GeV) 

Figure 2. The running gluon condensate as a function of Q when extracted from the 
Adler function. The grey area represents the allowed region using perturbation theory 
at 0(af) and the striped region using perturbation theory at ©(a,). The thick solid 
line denotes the running gluon condensate (CI = Q) from the lattice-inspired ansatz. 



for asymptotically large Q [ 13]. We have checked that the known perturba- 
tive contributions to the Adlcr function show good convergence properties. 
For Q between 1.0 and 2.0 GeV and setting ^ — Q, the 3-loop (two-loop) 
corrections amount to 5% (7%) and 0.5% (1.4%), respectively. On the other 
hand, the 0{as) correction to the Wilson coefficient of the gluon condensate 
are between —16% and —6% for /i between 0.8 and 2.0 GeV. Wc have com- 
pared the local dimension 6 contributions contained in the running gluon 
condensate with the sum of all dimension 6 terms in the OPE as they were 
determined in [ 13]. We again found that the corresponding dimension 6 
contributions have equal sign and roughly the same size. 

Our result for the running gluon condensate as a function of Q is shown 
in Fig. 2. for 1 GeV< Q < 2 GeV. The uncertainties are due to the ex- 
perimental errors in the spectral function and a variation of the renormal- 
ization scale /j. in the range Q ± 0.25 GeV. The strong coupling has been 
fixed at as{Mz) = 0.118. We have restricted our analysis to the range 
1 GeV < Q < 2 GeV because for Q < 1 GeV perturbation theory be- 
comes unreliable and for Q 2 GeV the experimentally unknown part 
of the spectral function at s > 2.2 GeV^ is being probed. The thick black 
line in Fig. 2 shows the running gluon condensate obtained from the lattice- 
inspired model of the last section for O = Q. We find that our result for the 



9 



running gluon condensate as a function of Q is consistent with a function 
that increases with Q. 

3.2.2. Charmonium sum rules 

The determination of the local gluon condensate G^^(oo) from charmo- 
nium sum rules was pioneered in Refs. [15]. By now there is a vast literature 
on computations of Wilson coefficients to various loop-orders and various 
dimensions of the power corrections in the correlator of two heavy quark 
currents including updated analyses of the corresponding sum rules (see 
e.g. Rcfs. [ 16,17]). 

The relevant correlator is 

where jfj, = cry^c^ and the n-th moment is defined as 

•^"nA- w) "°<«^' ■ 

Assiiming analyticity of 11'^ in away from the negative, real axis and 
employing the optical theorem, the nth moment can be expressed as a 
dispersion integral over the charm pair cross section in e+e~ annihilation, 

where s is the square of the cm. energy and Qc = 2/3. We consider the 
ratio [ 15] 

and extract the running gluon condensate as a function of n from the equal- 
ity of the theoretical ratio using Eq. (19) and the ratio based on Eq. (20) 
determined from experimental data. 

For the experimental moments we use the compilation presented in Ref. [ 
18], where the spectral function is split into contributions from the char- 
monium resonances, the charm threshold region, and the continuum. For 
the latter the authors of Ref. [ 18] used perturbation theory since no ex- 
perimental data are available for the continuum region. We have assigned 
a 10% error for the spectral function in the continuum region. For the 
purely perturbative contribution of the theoretical moments we used the 
compilation of analytic 0{a1) results from Ref. [ 18] and adopted the MS 
mass definition mc{rnc) (for any renormalization scale n). We have used 
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Figure 3. The running gluon condensate as a function of n when extracted from the 
ratio of moments Tn for rrlcijnc) = 1.23 (white triangles), 1.24 (black stars), 1.25 (white 
squares) and 1.26 GeV (black triangles). The area between the upper and lower symbols 
represent the uncertainties. The thick solid line indicates the running gluon condensate 
as it is obtained from the lattice-inspired model for g{x).) 

the one-loop expression for the Wilson coefficient of the gluon condensate, 
and wc have checked that for n < 8 the perturbative 0(q!^) corrections do 
not exceed 50% of the 0(q;.,) corrections for ji between 1 and 4 GeV. 

Our result for the running gluon condensate as a function of n is shown 
in Fig. 3. For each value of mc(mc) and n the area between the upper and 
lower symbols represents the uncertainty due to the experimental errors of 
the moments and variations of as{Mz) between 0.116 and 0.120 and of ji 
between 1 and 3 GeV. We see that the running gluon condensate appears to 
be a decreasing function of n. Since the exact width of the short-distance 
function / for the moments is unknown we estimate it using physical ar- 
guments. For large n, i.e. in the nonrelativistic regime, the width is of the 
order of the quark cm. kinetic energy mv^, which scales like m/n because 
the average quark velocity in the n-th moment scales like 1/^/n. [ 19] For 
small n the relevant short-distance scale is just the quark mass. Therefore, 
we take Vt = 2mc/n as the most appropriate choice of the resolution scale. 
In Fig. 3 the running gluon condensate is displayed for — 0.7 GeV 
and, exemplarily, Aq = 0.04 (i.e. {{as/T:)G'^)ut{oo) = 0.024 GeV*) and for 
n = (2.5 GeV)/n as the thick solid hue. We find that the n-dependence 
of the lattice inspired running gluon condensate and our fit result from the 
charmonium sum rules are consistent. However, the uncertainties of our ex- 
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traction are still quite large, particularly for n = 2, where the sensitivity to 
the error in the continuum region is enhanced. We also note that the values 
for the running gluon condensate tend to decrease with increasing charm 
quark mass. For mc{mc) = 1.27 GcV wc find that {(as/TT)G'^) ranges from 
—0.01 GeV^ to +0.01 GeV^ and shifts further towards negative values for 
rficirnc) > 1.28 GeV for all n > 2. Assuming the reliability of the DOE 
as well as the local OPE in describing the nonpcrturbativc effects in the 
charmonium sum rules and that the Euclidean scalar functions D[x) and 
Di{x) ill the parametrization (14) are positive definite, our result disfavors 
mc(mc) Z 1-28 GeV. 

3.3. A model calculation 

Among the early applications of the OPE in QCD was the analysis of 
nonperturbative effects in heavy quarkonium systems [ 20,21,22]. Heavy 

quarkonium systems arc nonrelativistic quark-antiquark bound states for 
which there is the following hierarchy of the relevant physical scales m 
(heavy quark mass), mv (relative momentum), mv^ (kinetic energy) and 
Aqcd: 

m > mt; > mv^ > Aqcd • (22) 

Thus the spatial size of the quarkonium system ~ {mv)~^ is much smaller 
than the typical dynamical time scale ^ {mv'^)^^ effectively rendering the 
problem 1-dimensional. In this section we demonstrate the DOE for the 
nonperturbative corrections to the n^^+^Lj = 1^5i ground state. We 
adopt the local version of the multipole expansion (OPE) for the expansion 
in the ratios of the scales m, mv and mv"^. The resolution dependent 
expansion (DOE) is applied with respect to the ratio of the scales mv^ 
and Aqcd- The former expansion amounts to the usual treatment of the 
dominant perturbative dynamics by means of a nonrelativistic two-body 
Schrodinger equation. The interaction with the nonperturbative vacuum is 
accounted for by two insertions of the local xE dipole operator, E being 
the chromoelectric field [ 20] . The chain of VEV's of the two gluon operator 
with increasing numbers of covariant derivatives times powers of quark- 
antiquark octet propagators [ 20] , i.e. the expansion in A/mv^, is treated 
using the DOE. 

At leading order in the local multipole expansion with respect to the 
scales m, mv, and mv"^ the expression for the nonperturbative corrections 
to the ground state energy reads 

/oo 
dtf{t)g{t), (23) 
-OO 



12 



where 



(24) 



with 



a,Uv.-i-] = ^(2fc,T)'(2%)'e-»('+») 

\ TO / ZTT 

L2'+i(2fc.T) L2(+i(2fcy)s! 



E 



(s + ; + i-^^)(s + 2/ + i)!' 



s=0 ' ' 12 fe 

?!>(a;) = ^ e"^^ , A; = - too^ . (25) 

The term Go is the quark- antiquark octet Green- function [ 22] , and cj) de- 
notes the ground state wave function. The functions P„ and are the 
Legendre and Laguerre polynomials, respectively. We note that t is the 

Euclidean time. This is the origin of the term it in the exponent appearing 
in the definition of the function /. In Fig. 4 the function fit) (solid line) 
is displayed for to = 5 GeV and ag = 0.39. Since the spatial extension of 
the quarkonium system is neglected and the average time between inter- 
actions with the vacuum is of the order of the inverse kinetic energy, the 
characteristic width of / is of order {mv'^)~^ ~ {ma1)~^ . As a comparison 
we have also displayed the function [J dt' f{t')eQ{t')]e2{t) for f2 = k'^/m 
(dashed line), which is the leading term in the delocalized multipole ex- 
pansion of /. The values of the first few local multipole moments /„(oo), 
which correspond to the local Wilson coefficients, read 

/o(oo) = 1.6518 3^1^, /2(^) = 1.3130 gf^, 

/4(oo) = 7.7570 3f^, /6(oo) = 30.4923!^, (26) 

/8(oo) = 4474.1 , /lo(oo) = 262709.3 ^ 

The term /o(oo) agrees with Ref. [ 21,22] and /2(oo) with Ref. [ 14]. The 
results for /ri>2(oo) are new. 

Let us compare the local expansion of E"p with the resolution- 
dependent expansion using the basis functions of Eq. (9). For the nonper- 
turbative gluonic field strength correlator we use a lattice inspired model 
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Figure 4. The perturbative short-distance function f{t) for m = 5 GeV and as = 0.39 
(soUd hne) and the leading term in the delocaUzed multipole expansion of f(t) for U = 
k^/m (dashed hne). 

of the form 

git) = 12 Ao exp ( - ^J^^JXa + l) , 

^0 = 0.04 GcV'' , A;^^ ^ 0.7 GeV , (27) 

This model has an exponential large-time behavior and a smooth behavior 
for small t. The local dimension 4 gluon condensate in this model is 

G»,G^,) = ^ = 0.024 GeV^ . (28) 

In Tab. 1 we have displayed the exact result and the first four terms of 
the resolution-dependent expansion of E^^ for the quark masses m = 
5, 25, 45, 90, 175 GeV and for Q = oo and O = k^/m. For each value of the 
quark mass the strong coupling has been fixed by the relation = as{k). 
Note that the scries arc all asymptotic, i.e. there is no convergence for any 
resolution. The local expansion (Q = oo) is quite badly behaved for smaller 
quark masses because for fc^/m ^ A^^ any local expansion is meaningless. 
In particular, for m = 5 GcV the subleading dimension 6 term is already 
larger than the parametrically leading dimension 4 term. This is consis- 
tent with the size of the dimension 6 term based on a phenomenological 
estimate of the local dimension 6 condensate. For quark masses, where 
fc^/m ^ A^^, the local expansion is reasonably good. However, for finite 
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Table 1. Nonperturbative corrections to the heavy quarkonium ground state level 
at leading order in the multipole expansion with respect to the scales m, mv 
and mv^ for various quark masses m based on the model in Eq. (27). Dis- 
played are the exact result and the first few orders in the DOE for £7 = oo and 

n = k-hii. The iiunilKTs an- rounded oil In uuils of 0.1. 0.01 or 0,001 MoV. 





Q = CO 


n = k'^/m 


(GcV) 


OLa 


j TTl 

(McV) 


(McV) 


n 


fn9n 

(MeV) 


y^" fa 

^i— j^y^ 

(MeV) 


fnQn 

(McV) 


fa 

^i— J^y^ 

(McV) 


5 


0.39 


0.338 


24.8 





38.6 


38.6 


24.2 


24.2 


2 


-65.7 


-27.2 


-3.9 


20.3 


4 


832.7 


805.5 


12.1 


32.4 


6 


-35048.0 


-34242.4 


-43.1 


-10.8 


25 


0.23 


0.588 


12.6 





16.0 


16.0 


12.8 


12.8 


2 


-9.0 


7.0 


-1.2 


11.5 


4 


37.8 


44.8 


2.6 


14.1 


6 


-526.6 


-481.8 


-6.7 


7.4 


45 


0.19 


0.722 


4.9 





5.9 


5.9 


5.0 


5.0 


2 


-2.2 


3.7 


-0.4 


4.6 


4 


6.1 


9.8 


0.7 


5.3 


6 


-56.4 


-46.6 


-1.4 


3.8 


90 


0.17 


1.156 


1.05 





1.15 


1.15 


1.07 


1.07 


2 


-0.17 


0.98 


-0.04 


1.02 


4 


0.18 


1.16 


0.04 


1.07 


6 


-0.65 


0.51 


-0.07 


1.00 


175 


0.15 


1.750 


0.245 





0.258 


0.258 


0.249 


0.249 


2 


-0.016 


0.242 


-0.005 


0.244 


4 


0.008 


0.250 


0.003 


0.247 


6 


-0.012 


0.237 


-0.003 


0.244 



resolution Q. = jm, the size of higher order corrections is considerably 
smaller than in the local expansion for all quark masses, and the series 
appears to be much better behaved. The size of the order n term is sup- 
pressed by approximately a factor 2~" as compared to the order n term in 
the local expansion. We find explicitly that terms in the series with larger 
n decrease more quickly for finite resolution scale as compared to the local 
expansion. One also observes that even in the case fc^/m < A^^, where 
the leading term of the local expansion overestimates the exact result, the 
leading term in the delocalized expansion for O = fc^/m agrees with the 
exact result within a few percent. It is intuitively clear that this feature is 
a general property of the delocalized expansion. 

We would like to emphasize the above calculation is not intended to 
provide a phenomenological determination of nonperturbative corrections 
to the heavy quarkonium ground state energy level, but rather to demon- 
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strate the DOE within a specific model. For a reahstic treatment of the 
nonperturbative contributions in the heavy quarkonium spectrum a model- 
indcpcndcnt analysis should be carried out. In addition, higher orders in 
the local multipole expansion with respect to the ratios of scales m, mv and 
mv^ should be taken into account, which have been neglected here. These 
corrections might be substantial, in particular for smaller quark masses. 
However, having in mind an application to the bottomonium spectrum, we 
believe that our results for the expansion in A/mv^ indicate that going 
beyond the leading term in the OPE for the bottomonium ground state 
is probably meaningless and that calculations based on the DOE with a 
suitable choice of resolution are more reliable. 

4. Summary and Outlook 

In this talk I have presented a generalization of the expansion of the cor- 
relator of gauge invariant currents into local operator averages. It is built 
on appropriate choices of basis sets spanning the dual spaces which corre- 
spond to the interplay between perturbative and nonperturbative Af-point 
correlations in the Euclidean propagation of a gauge invariant current. It 
was demonstrated that this framework, when applied to the V + A and 
charmonium systems, yields an experimentally confirmed running gluon 
condensate. In a model calculation for the nonperturbative corrections to 
the ground state of a heavy quarkonium superior convergence properties of 
the DOE were obtained when compared with the OPE. 

The framework has a wealth of potential applications. For example, it 
can be used to improve the scaling relation for decay constants of heavy- 
light mesons as it is obtained in HQET (see [ 9] for an attempt). The issue of 
the violation of local quark-hadron duality by the OPE can be re-addressed 
[ 23], and a scheme can be developed to estimate duality violation effects. 
This is of particular importance for the determination of standard model 
parameters from mixing and decay of B mesons [ 6] . 
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